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FOREWORD 


As  part  of  the  on-going  research  program  in  "Decision  Control 
Models  in  Operations  Research,"  Mr.  J.  Christopher  Mitchell  examines 
the  problem  of  minimizing  costs  in  a  multi-item  system  by  an  appro¬ 
priate  choice  of  single-item  service  objectives.  Mr.  Mitchell  inves¬ 
tigates  the  properties  of  optimal  policies  and,  in  the  special  case 
of  exponentially  distributed  demand,  develops  an  efficient  algorithm 
for  computing  optimal  policies.  He  shows  that  this  approach  can 
lead  to  significant  savings  over  the  uniform  service  approach  that 
is  often  used  in  applied  settings. 

This  report  is  somewhat  of  a  departure  from  earlier  reports  in 
that  it  is  the  first  to  consider  a  system  of  interdependent  items. 
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1.  INTRODUCTION 


We  consider  a  simple  multi-item  inventory  system  and  study  the 
problem  of  specifying  system-wide  service  objectives.  Most  theoret 
i cal  models  considered  in  the  literature  assume  that  holding  or 
shortage  costs  are  applied  to  any  excess  inventory  or  unsatisfied 
demand,  respectively.  A  major  difficulty  in  applying  such  models 
in  practice  is  in  the  specification  of  shortage  costs.  Frequently, 
the  manager  will  set  shortage  cost  parameters  based  on  an  objective 
of  satisfying  demand  with  at  least  some  minimum  probability.  Such 
an  approach  is  often  preferred  because  subjective  factors  can  more 
easily  be  expressed  as  a  probability  of  satisfying  demand  than  as 
a  cost  for  each  shortage  incurred. 

A  significant  shortcoming  of  this  approach  is  that  it  usually 
entails  the  setting  of  one  system-wide  probability  of  demand  satis¬ 
faction  that  is  then  applied  uniformly  to  all  items  in  the  system. 
In  this  paper  we  raise  the  issue  of  specifying  different  service 
objectives  for  individual  items  while  still  satisfying  some  given 
system-wide  objective.  The  value  of  this  approach  is  that  system 
costs  can  be  reduced  below  those  of  a  method  that  requires  identi¬ 
cal  service  for  all  items. 

Specifically,  we  consider  a  base-stock  inventory  model,  a  mult 
item,  multi-period  inventory  problem  with  stochastic  demands,  fixed 
leadtimes,  and  stationary  ordering,  holding,  and  shortage  costs. 

For  each  item  i  there  is  an  ordering  cost  c^  per  unit  ordered, 
incurred  upon  delivery  L.  2  0  periods  after  order  placement. 

There  Is  also  a  holding  cost  h^  per  unit  of  inventory  on  hand, 
incurred  at  the  end  of  each  period.  The  demand  realized  for  item 
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i  in  period  n  is  Din,  with  absolutely  continuous  cdf  4^  and 
density  4^.  We  assume  that  Din  2  0  with  probability  1*  and 
♦.  (0)  =  0.  The  mean  of  4^  is  y^,  and  all  demands  D^n  are 
mutually  independent.  The  initial  inventory  level  for  item  i  is 
ui ,  and  the  single  period  discount  factor  is  B,  where  0  <  B  <  1. 
There  are  N  items  and  T  (possibly  infinite)  periods  in  which 
decisions  must  be  made.  When  T  is  finite,  there  are  also  final 
salvage  periods  T  +  +  1  for  i  =  1,...,N. 

A  solution  to  the  single-product  version  of  this  problem  is 
given  by  Veinott  [2,  pp.  754-757]  when  a  unit  shortage  cost  pi , 

incurred  at  the  end  of  each  period,  is  specified  for  each  item  i. 

★ 

Let  denote  the  (L^U-fold  convolution  of  4>.j  with  itself. 

There  is  a  myopic  optimal  policy  x  =  (x-j . x^)  given  by 

*  -  ( 1  - 6 ) c i 

= - -  - - ,  where  is  the  base-stock  level  for  item 

pi  ni 

i  in  all  periods.  Thus  the  policy  requires  that  whenever  the  in¬ 
ventory  position  (stock  on  hand  plus  stock  on  order)  for  item  i 
falls  below  x^,  an  order  is  placed  to  increase  it  up  to  x^ .  The 
probability  of  meeting  the  demand  for  any  period  for  item  i  is 
Pi  -  (l-S)c. 

thus  — — - — ~u - .  Typically,  a  manager  wishes  to  meet  all  demands 

pi  1 

in  all  periods  with  a  probability  of  at  least  a,  for  some  pre¬ 
scribed  0  <  a  <  1,  and  hence  sets  each  shortage  cost  pi  such 


that 


Pi  -  (1-S)ci 
i>i  *  hi 


=  a. 


This  paper  will  discuss  an  alternative  base-stock  inventory 
model,  one  without  shortage  costs.  Instead  of  these  costs,  we 
require  that  some  specified  set  of  service-level  constraints  be 
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satisfied  in  each  period.  Except  for  the  most  simple  constraints, 
it  is  considerably  more  difficult  to  find  an  optimal  policy  when 
service-level  constraints,  rather  than  shortage  costs,  are  used. 

We  formulate  this  model  as  a  constrained  nonlinear  program 
(NLP).  We  discuss  two  types  of  service-level  constraints.  The 
first  is  the  one  frequently  used  by  managers,  namely,  requiring 
that  all  demands  in  all  periods  be  met  with  a  probability  of  at 
least  a.  The  second  is  more  general  than  the  first.  It  requires 
that  all  demands  for  item  i  in  all  periods  be  met  with  a  proba¬ 
bility  of  at  least  a. ,  for  i  =  1,...,N,  where  the  ai  are  not 
necessarily  all  the  same.  We  also  require  a  system-wide  constraint 
on  the  ai  in  each  period.  Each  item  is  controlled  by  a  base-stock 
policy  that  assures  demand  satisfaction  with  probability  for 

item  i.  Moreover,  a  solution  x  =  (x^,...,xN)  of  this  constrained 
NLP  is  equivalent  to  specifying  shortage  costs  pi  for  each  item 
p.  -  (l-B)c. 

i  so  that  - ^-rj - =  a..  Although  we  assume  that  the  demand 

Pi  .K  i 

distributions  are  absolutely  continuous,  we  note  that  if  the  demand 
distributions  are  discrete  a  deterministic  policy  may  not  be  opti¬ 
mal  [2,  p.  757].  Nonetheless,  we  would  recommend  that  research  be 
focused  on  deterministic  base-stock  policies,  since  these  are  more 
appropriate  in  practice. 

For  our  constrained  NLP  we  derive  equations  representing  first- 
order  necessary  and  second-order  sufficient  conditions,  the  solution 
to  which  is  an  optimal  policy.  However,  it  is  not  clear  how  to 
solve  these  equations  efficiently  for  a  general  demand  distribution. 
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For  exponential  demand  distributions  we  derive  a  computationally 
efficient  algorithm  to  solve  these  equations.  We  then  consider  an 
example  in  which  we  compare  the  total  system  costs  when  each, of  the 
two  types  of  service-level  constraints  is  used.  We  show  that  there 
are  inventory  systems  of  practical  interest  for  which  there  is  a 
substantial  cost  savings  when  using  the  general  constraint  rather 
than  the  simpler  one  often  used  by  managers. 

2.  THE  GENERAL  INVENTORY  MODEL 

Recall  that  for  our  alternative  base-stock  model  we  have  for 
each  item  i  =  1,...,N  and  each  period  n  =  1,...,T+L.,  the  fol¬ 
lowing:  ordering  cost  ci ,  holding  cost  hi ,  fixed  leadtime 
L.  >  0,  demands  0.  with  continuous  cdf  <f . ,  density  41.  and 
mean  ,  initial  inventory  ,  and  a  single-period  discount  fac¬ 
tor  6  with  0  <  8  <  1.  Each  item  i  is  controlled  by  a ‘myopic 
base-stock  policy  x^ ,  and  x  =  (x^,...,xN).  We  seek  to  minimize 
the  total  expected  discounted  cost  over  the  entire  horizon  subject  to 
a  system-wide  service  constraint  being  satisfied  in  every  period. 

Let  Cni(x.j)  be  the  expected  discounted  cost  for  item  i  in 
period  n  +  Li .  If  the  leadtime  is  greater  than  0,  costs  in 
periods  l,...,Li  will  be  independent  of  x^ .  Now  for  n  =  1,...,T 
and  i  =  1 , . . . ,N, 

n+L.-lp  1 

Cni<»i>se  ♦'iE<Vj  •  <» 

where  Y  .  is  the  inventory  on  hand  for  item  i  at  the  end  of 
mi 

period  m,  and  Rn-  is  the  order  size  for  item  i  placed  in 
period  n  (and  hence  received  and  paid  for  in  period  n  +  L.). 
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Since  is  inventory  on  hand  plus  on  order  in  each  period, 

/  n+L,  \+ 


where  (y)+  =  max^O.y}.  And  since  the  policy  is  a  base-stock 
pol i cy , 


Rin  = 


xi  -  ui  n  =  1 


n  =  2.....T 


Note  that  the  value  of  1  implies  that  if  the  initial  inventory 
u.j  exceeds  the  base-stock  level  xi ,  the  excess  is  salvaged  at  the 

unit  ordering  cost  ci . 

★ 

Recall  that  ♦  .  is  the  (L^+lJ-fold  convolution  of  *.  with 
★ 

itself.  Let  4>.  be  its  density.  Then  since  the  D . .  are  iid 
with  cdf  'I>i  for  j  =  n,...,n+L.j»  we  have 


n+L.  v+ 


('  Li  \ 


xi 

=  ^  (xi-t)d«>i (t)  . 


Since  ^(0)  =  0  for  all  i,  integration  by  parts  yields 


fi  * 

E(Yn+L  fi)  =  J  4>i(t)dt  for  i  =  1 . N  , 


n  =  1,...,T  . 


We  also  have 
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E<Rni>  = 


x.  -  u, 


(3) 


We  assume  a  salvage  operation  for  item  i  in  period 
T  +  L-  +  1,  where  any  excess  stock  is  sold  at  its  unit  ordering 
cost,  and  any  unfilled  demand  is  satisfied  immediately.  Hence, 
the  salvage  cost  is  given  by 


which  has  expected  value 

T+Li 

CT+1 , i ^ xi ^  =  ~6  ci^xi  ‘  (Li+1)ui]  • 

Our  introduction  of  salvage  cost  follows  the  development  in  [2], 
where  it  is  on  artifice  to  guarantee  a  myopic  optimal  policy. 

When  the  horizon  T  is  long,  the  salvage  terms  are  an  acceptable 
approximation.  For  the  infinite  horizon  case,  which  is  considered 
at  the  end  of  Section  3,  the  artifice  evidently  vanishes. 

T+l  N 

let  C(x)  =  7  T  C  - (x . ) ,  the  total  discounted  expected 

n=l  i=l  n1  1 

cost  excluding  the  policy-independent  costs  incurred  in  periods 
1,...,L,  where  L  =  min{L^:i  =  1,...,N).  Also  let 


H(x)  = 


n  .  x. 


(t)dt 


and 


L. 


h(x)  e  1c.xi 

i*1  1  1 


Then  (1),  (2),  (3),  and  (4)  imply  that 


T  N  n+|  ,  x.i 

C(X)  .  E  E  3  r  h 

n=l  i=l  0 


i  ^  4>*(t)dt 


N 

+  E 

i=l 


L-  E  n+L.-l 

3  c^x.-u.)  +  Lj  3  c-y. 
111  n=2  11 


T+L. 

-  E  3  ci[xi  -  (L.+l)vi]  . 


After  some  algebraic  simplifications,  for  6  <  1  we  can  write 


C(x)  =  C(x)  +  (^)E  B^c^  +  6T  E  eLic.(l1.+l)yi 


i*l 


(5) 


where  C(x)  =  ij~-{H(x)  +  (l-g)h(x)]  .  (6) 

Note  that  C(x)  -  C(x)  are  policy-independent  costs. 

In  sunwary,  we  have  that  the  total  expected  discounted  cost 
excluding  all  policy- independent  costs  is 


C(x)  =  y^-[H(x)  +  c(x)]  , 


(7a) 


gyp  *fwq 
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(7b) 

(7c) 

We  are  interested  in  minimizing  C(x)  subject  to  a  service 
constraint  on  the  probability  of  satisfying  demand.  We  consider 
two  types  of  service  constraints.  The  first  type  requires  that 
the  demand  for  every  item  in  every  period  be  satisfied  with  a  prob¬ 
ability  of  at  least  a,  where  0  <  a  <  1.  We  call  this  service 
requirement  Independent  and  Identical  Service  (IIS).  Thus,  there 
will  be  N  service  constraints 

4>*(x.)  *  a  for  i  *  1,.. .  ,N  ,  (8) 

which  must  be  satisfied  in  every  period.  IIS  is  an  approach  fre¬ 
quently  used  in  practice.  Indeed,  if  we  make  the  reasonable  as¬ 
sumption  that  for  some  myopic  optimal  policy  x  =  (x, ,...,x^) 
these  constraints  are  tight,  finding  such  a  policy  is  equivalent 
to  specifying  the  shortage  costs  p^  in  Veinott's  model  [2]  des¬ 
cribed  in  the  Introduction  by 

Pi  -  (l-B)ci 

— p  — =  a,  for  i  =  1,...,N  . 

The  second  type  of  service  constraint  Is  a  generalization  of 
IIS  and  is  called  General  Service  (GS ) .  For  this  service  the  sys¬ 
tem  of  N  constraints  on  the  probability  of  satisfying  demand  is 


where  H(x) 


N  .  L  .  xi 

*  \))  ♦;<*>« 


N  L 

and  c{x)  =  (l-S)h(x)  =  2  (1-6) (6  ’cjx* 

i=l 


i  i 
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replaced  by  a  single  constraint  requiring  that  a  weighted  average 
of  the  probabilities  be  at  least  o.  The  value  of  this  approach 
is  that  varying  service  levels  can  be  applied  to  the  items  se  as 
to  reap  the  greatest  benefits  in  system-wide  costs.  We  use  the 
constraint 


N 

2  •  o]  i  o  , 

i=l  1  1  1 


(9) 


where  k^,...,kN  >  0  are  specified  weights.  In  the  following 
sections  we  often  set  k^  =  y^. 

We  will  require  that  constraint  (9)  be  satisfied  in  every 
period.  This  does  not  destroy  the  myopic  character  of  the  optimal 
solution. 

We  minimize  the  cost  function  given  by  (7a),  (7b),  and  (7c) 
subject  to  the  constraint  (9).  In  summary,  we  have  formulated  the 
following  optimization  problem: 


minimize  C(x)  =  ^ +  C(*)J 


subject  to  g(x)  *  X*  k.[i>. (x.)  -  a]  z  0  , 
i=l  1  1  1 


N  L  * '  i 

where  H(x)  =  S  8  'h.  \  4>.(t)dt 

i=l  1  i  1 


L 

and  c(x)  =  X3  (1-6)6  V.,x. 

i=l  1  1 


(10) 


10 


We  mention  again  that  a  myopic  optimal  policy  x  =  (x^,...,xN) 
found  to  be  the  solution  of  this  NLP  is  equivalent  to  specifying 

P-  -  0-B)c.  * 

shortage  costs  p.  as  — —  —  =  ^(x^)  in  Veinott's  model 

[2]  discussed  in  the  Introduction.  Also,  the  probability  ai 
that  the  demand  for  item  i  be  met  in  any  period  is 
P^  -  0  - 6 ) c i 


3.  ANALYSIS 


3.1  First-  and  Second-Order  Optimality  Conditions 

We  use  NLP  techniques  to  analyze  (10).  Specifically,  we  de¬ 
rive  first-order  necessary  and  second-order  sufficient  conditions 
for  x  to  be  a  solution  of  (10).  We  also  investigate  when  this 
local  optimum  can  be  known  to  be  a  global  optimum. 

It  is  clear  that  if  (10)  has  a  solution  x,  it  has  one  with 
g(x)  =  0,  and  so  we  assume  that  g(x)  =  0  at  optimality.  Since 
there  is  only  one  constraint  g(x)  ;>  0,  and  g  is  differentiable 
everywhere  on  (0,+®),  a  first-order  necessary  condition  for  x 
to  solve  (10)  is  that  there  exist  A  2  0  such  that 


VC(x)  -  XVg(x)  =  0  , 
where  g(x)  =  0  , 


(ID 


as  proved  in  [1,  p.  41].  Now 


C(x) 


0 


fV[s  1hiVxi 


)  +  (1-6)6 
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and 


j|-9(x)  -  k^tx,)  • 
Hence  (11)  is  equivalent  to 


★ 

X 


/i-e\x  . 

''-67  (*i ) 


being  independent  of  i  for  all  i  =  1,...,N. 

Suppose  that  at  optimality  <f>*  and 
following  form  for  i  =  1,...,N: 

★ 

^(x^)  =  v^a,  sane  0  <  <  1/a 


.  el<[Vi»*  d-8)cj] 


A*k, 


some 


can  be  written  in  the 


X*  >  0  . 


(12) 


Since  X*  is  independent  of  i,  (12)  together  with  g(x)  =  0, 
which  is 


N 

E  k.[v,  -  1]  =  0 

i*l  1  1 


satisfy  the  necessary  condition  (11).  We  summarize  these  results  in 
the  following  theorem. 

Theorem  1:  A  first-order  necessary  condition  for  x  to  solve 
(10)  is  that  there  exist  v^,...,vN, 

X*  >  0,  such  that  for  i  =  1,...,N, 


with  0  <  v^  <  1/a,  and 


(a)  4>*(x.) 

(b)  *;(xf) 

N 

(C)  E  k^-lj 
i=l  1  1 

Given  that  the  first-order  condition  (11)  holds,  we  proceed  to 
derive  a  second-order  sufficient  condition  that  x  solve  (10).  As 
proved  in  [1,  p.  48],  a  second-order  condition  is  the  following: 

For  all  z  e  with  z  f  0, 

zT[V2C(x)  -  XV2g(x)]z  >  0  ,  (14) 

where  A  >  0  is  the  same  parameter  as  given  in  (11). 

Now 

_L1  * 

8  hj<t>j(xj)  0 

V2C(x)  =  *•  ^  ,  (15) 

0  6  ^  XN )- 

and 

(x^)  0 

729(x)  =  ^  .  (16) 

-  0  kN*N'(xN> 
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and  so  (14)  is  equivalent  to 


>  o  , 


1  ft 

because  X  =  X*. 

Since  it  is  sufficient  that  each  term  in  the  summand  be 
positive,  we  have  established  the  following  theorem: 


Theorem  2.  Given  that  (11)  holds  for  x,  a  second-order 
sufficient  condition  for  x  to  solve  (10)  is  that  for  i  =  1 . N, 


dXj  “MV  *  6  1  x*kn.  W  • 


(17) 


Note  that  (17)  is  trivially  satisfied  if  <J>^  is  positive  and 


nonincreasing  at  x^. 


3.2  Undiscounted  and  Infinite  Horizon  Models 

We  briefly  consider  what  changes  are  necessary  to  model  the 
systems  that  are  undiscounted  and/or  have  an  infinite  planning 
horizon.  It  turns  out  that  although  the  formulation  (10)  must  be 
modified  somewhat,  the  first-  and  second-order  conditions  (13)  and 
(17)  remain  intact. 

The  undiscounted  finite  horizon  model  is  obtained  by  setting 
8*1.  The  only  change  necessary  in  the  formulation  occurs  with 
respect  to  (7a),  where  it  is  assumed  B  <  1.  However,  when  8=1, 
(7a)  becomes 


C(x)  =  T[H(x)  +  c ( x ) ]  . 


Moreover,  when  8=1,  (7b)  and  (7c)  become 


H(x)  = 


and  c(x)  =  0,  respectively. 

Therefore  the  undiscounted  finite  horizon  model  analogous  to 
(10)  is 

\ 


minimize  C(x)  =  TH(x) 
N 


subject  to  g(x)  =  k.[$.(x.)  -  a]  >  0  , 

i*l  1  1 


N  *i  # 

where  H(x)  =  h.\i,(t)dt  . 

i=i  i 


(18) 


It  is  straightforward  to  verify  that  the  first-order  necessary 
and  second-order  sufficient  conditions  for  optimality  are- still 
(13)  and  (17),  respectively,  with  8  =  1  inserted  in  those  formulas. 

The  undiscounted  infinite  horizon  model  is  obtained  by  setting 
8  =  1  and  T  =  +®.  It  is  clear  from  (18)  that  the  total  expected 
undiscounted  cost  is  infinite,  and  so  we  must  use  an  alternative 
optimality  criterion.  We  seek  to  minimize  the  average  expected  cost 
per  period.  Thus  this  is  the  previous  model  (18)  with  T  =  1,  and 
so  the  optimality  conditions  (13)  and  (17)  hold  with  8  =  1  and 
T  *  1  (not  T  *  +*>)  inserted  in  those  formulas. 

The  discounted  infinite  horizon  model  is  obtained  by  setting 
T  =  -hc  (and  8  <  1).  In  the  formulation  of  the  model  (10)  the 
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only  dependence  on  T  is  found  in  the  objective  function,  which 
becomes 


C(x)  =  T~g£H(x)  +  c(x)3 

when  T  *  +».  In  order  to  guarantee  that  the  total  expected  dis¬ 
counted  cost  is  finite,  we  must  verify  that  the  policy-independent 
terms  in  (5)  are  finite  for  T  =  +®.  The  only  dependence  on  T 
are  the  terms 

TT  a"d  sT  • 


and  for  6  <  1 , 


lim  =  yto  <  and  lim  BT  =  0  . 

T-KO  C  C  J-¥00 

Thus  the  total  expected  discounted  cost  is  finite,  and  the  analog 
to  the  model  (10)  for  the  discounted  infinite  horizon  model  is 

minimize  C(x)  = 


subject  to  g(x)  =  52  k.[$.(x.)  -  a]  s  0  , 

1-1  1 


where  H(x) 


N  L  Xfi 

'  ”  I  A 


(t)dt 


and  c(x)  =  52  (1-6)8  ^c.x. 

1-1  1  1 


(19) 
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Again,  the  first-  and  second-order  optimality  conditions  (13) 
and  (17)  remain  unchanged,  as  is  easily  shown. 

3.3  Global  Optima 

We  now  seek  to  establish  conditions  that  guarantee  that  a 
local  optimum  to  (10)  is  in  fact  a  global  optimum.  In  general,  a 
local  optimum  need  not  be  a  global  optimum.  However,  if  (10)  is  a 
convex  NLP,  i.e.,  if  on  their  respective  domains  C(x)  is  a  convex 
function  and  g(x)  is  a  concave  function,  then  any  local  optimum 
is  also  a  global  optimum,  as  proved  in  [1,  p.  96].  Now  (15)  implies 
that  C(x)  is  always  convex  since  v  C(x)  is  always  positive  def¬ 
inite.  However,  (16)  implies  that  g(x)  is  concave  if  and  only  if 
★ 

is  nonincreasing  at  x^  for  all  i  =  1,...,N  (since  this  is 
2 

equivalent  to  v  g(x)  being  negative  semi-definite  at  x).  Recall 
that  this  is  almost  the  same  condition  that  guaranteed  that  the 
second-order  sufficient  condition  (17)  be  satisfied. 

In  our  examples,  we  consider  the  exponential  and  uniform  dis¬ 
tributions.  Since  their  density  functions  are  always  nonincreasing 
on  their  supports,  it  will  be  sufficient  to  search  for  local  optima 
using  (13)  and  (17),  for  they  will  also  be  global  optima. 

4.  EXPONENTIAL  DEMAND 

We  consider  the  special  case  where  Li  =  0  for  i  =  1,...,N, 
and  all  demands  are  exponentially  distributed.  Using  the  first- 
and  second-order  conditions  we  derive  a  computationally  efficient 
algorithm  to  compute  an  optimal  base-stock  policy  x. 


17 


Since  the  exponential  density  is  positive  and  decreasing  on 
its  support,  the  sufficient  condition  (17)  is  trivially  satisfied 
for  all  x  >  0.  Thus  we  need  only  search  for  x  >  0  which  satis¬ 
fies  the  necessary  condition  (13).  Because  =  0, 


^ (x)  =  ^ (x)  = 


1  .  e(‘x/Pi) 


x  <  0 

X  2;  0 


and 


i>*(x) 


<fhj(x)  = 


J_  gt-x/u^ 
yi 


x  <  0 

x  >  0 


So  when  x  >  0,  (a)  and  (b)  of  (13)  are  equivalent  to 


X*ki  -  ui(l-6)ci 
Vi  u(.\*ki  +  L^TT 


(20) 


This,  together  with  (c)  of  (13),  is  equivalent  to 


N 


E 

i=l 


ki 


X*ki  -  wi(l-B)ci. 
a  (X*ki~V^JThi") 


=  0  . 


Now  ^(x-j)  =  v^a  if  and  only  if 


(21) 


xi  =  -pi  log  (1  -  v^)  ,  (22) 

and  so  if  we  can  solve  (21)  for  X*  >  0,  (20)  gives  us  v..  and 

(22)  gives  us  x..  We  have  the  following  theorem: 
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Theorem  2:  There  is  a  unique  X*  >  0  that  solves  (21). 


Proof:  Let 


f(X)  = 


Ak.j  -  y,. (l-B)c^ 
ki  a(Xk^  +  y.h.)  1 


Then  f(0)  <  0  and  f{+°°)  >  0.  Since  f  is  continuous  on 
(0,+®),  the  Intermediate  Value  Theorem  implies  that  there  exists 
X*  c  ( 0 ,+o°)  such  that  f(X*)  =  0.  It  is  straightforward  to  verify 
that  f'(X)  >0  on  (0,+®),  and  so  Rolle’s  Theorem  implies  that 
X*  is  unique.  0 


Li 


^i: 


Thus  we  have  the  following  algorithm  to  compute  x  when 
0  and  the  demands  are  exponentially  distributed  with  means 


(1)  Solve 


Y"  .  pki  •  V'6)ci  .1  . 


for  X*  >  0  . 


X*k.  -  p.(l-6)Ci 

(2)  Set  for  1  ■ 1 . N- 


i  "li 


(3)  Set  x.  = 


.  =  -yi  log  (1-  v,.a)  for  i  =  1 .  ,N.  J 


(23) 


Step  (1)  can  be  done,  for  example,  using  a  bisection  method. 
Even  for  very  large  N,  this  is  not  computationally  very  difficult. 

Recall  that  in  the  Introduction  we  discussed  two  possible 
service  constraints  we  might  impose  on  the  inventory  system.  The 
first,  which  we  call  Independent  and  Identical  Service  (IIS), 
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requires  that  we  use  the  N  constraints  given  in  (8)  for  each 
period.  The  second  service  constraint,  which  we  call  General 
Service  (GS),  requires  that  we  use  the  single  constraint  (9)  in 
each  period. 

If  we  use  the  IIS  constraint,  the  fact  that  the  constraints 
are  tight  at  optimality  implies  that  the  optimal  value  for  x  can 
be  found  by  solving  the  system  of  equations 

*.j(x.)  =  a,  i  =  1 . N  , 

for  x  =  (x, . x^).  Hence,  the  optimal  base-stock  policy  is  given 

by 


xi  =  -yi  log(l-a)  .  (24) 

If  we  use  the  GS  constraint,  the  optimal  base-stock  policy  x 
must  be  found  using  the  algorithm  given  in  (23).  Hence,  the  impact 
of  imposing  the  GS  constraint  as  opposed  to  the  IIS  constraint  is 
to  replace  y.  log(l-o)  in  (24)  with  ui  log(l-v..a),  where  v^ 
is  algorithmically  computed. 

It  is  clear  that  the  GS  constraint  is  less  restrictive  than 
the  IIS  constraint,  and  hence,  the  optimal  objective  value  of  (7a) 
subject  to  the  GS  constraint  is  no  larger  than  that  subject  to  the 
IIS  constraint.  However,  it  is  a  good  deal  more  work  to  find  the 
base-stock  policy  in  the  GS  case  than  in  the  IIS  case.  One  may 
wonder  for  what  system  parameters,  if  any,  the  extra  work  and  ex¬ 
pense  involved  in  using  (23)  is  justified.  In  other  words,  when, 
if  ever,  is  the  cost  function  (7)  evaluated  at  the  solution  (23) 
of  (10)  significantly  smaller  than  when  evaluated  at  (24)?  We 
explore  this  question  with  an  illustrative  two-item  system. 
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We  consider  a  single  period,  two-item  system  that  is  undis¬ 
counted.  It  is  straightforward  to  verify  that  the  objective 
function  (10)  is 

2 

C(x)  =  H  h . [x .  -  y.<Mx.)3  •  (25) 

i=l  1  i  111 

We  denote  the  optimal  base-stock  policy  for  the  IIS  and  GS 
models  by  x  =  (xi .  ,x^)  and  x  =  (xlj1, . . .  ,x^) ,  respectively. 
Then  (24)  and  (23)  imply,  respectively,  that  for  i  =  1,...,N, 

xj  -  -yi  log  (l-o) 

(26) 

p 

xn-  =  log  (1-v.a) 

where  the  v.  are  found  by  algorithm  (23). 

Let  D  be  the  percent  decrease  in  minimum  system  cost  when  the 
GS  constraint,  rather  than  the  IIS  constraint,  is  used.  Then, 


D  =  D(y^  ,h^  ,  ) 


and  (25)  and  (26)  imply  that 
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0  =  100  1 


2 

[a  +  log(l-a)]  S  h.u. 

i=l  1  1 

_ 

Stv.a  +  log(l-v.a)]h.y. 
i=l  1  iii 


(27) 


In  order  to  investigate  the  sensitivity  of  D  to  the  param¬ 
eters  ui ,hi  and  a,  we  show  that  for  fixed  a,  D  is  a  function 
only  of  and  t^/h-j.  We  then  graph  D  evaluated  at  pairs 

(v^/ui  .f^/hi ) ,  and  see  for  what  pairs  D  is  large. 

Theorem  ?:  For  fixed  a,  D  as  given  in  (27)  is  a  function 
only  of  and  h^/h^ . 

Proof:  The  proof  is  in  the  Appendix.  0 

In  Figures  1  and  2  we  have  plotted  the  level  curves  of  D  for 
a  -  .8  and  a  =  .9,  respectively.  We  used  the  IMSL  routine 
ZBRENT  to  compute  Step  (1)  of  algorithm  (23). 

We  consider  0  =  5%  to  be  a  significant  decrease  in  optimal 
cost.  In  Figures  1  and  2  we  have  plotted  the  level  curves  corres¬ 
ponding  to  D  =  5.  We  call  the  two  regions  in  each  figure  for 
which  D  s  5  critical  regions.  Whenever  a  pair  (p2/Ti-|  »h2/,hl  ^ 
falls  within  a  critical  region,  a  significant  decrease  in  optimal 
cost  results  from  using  the  GS  rather  than  the  IIS  constraint.  The 
symmetry  observable  in  the  graphs  follows  from  the  easily  verifiable 
fact  that 
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Note  that  the  largest  values  of  D  are  found  in  the  regions 
where  f^/h^  and  are  large,  and  where  h^h^  and 

are  large.  A  system  for  which  this  is  typical  is  one  in  which  some 
items  have  a  very  low  expected  demand  and  a  very  high  holding  cost, 
or  vice  versa  [high  holding  cost  £  high  unit  purchase  cost 
(s  expensive  items)].  Many  inventory  systems  of  practical  interest 
contain  such  items.  For  such  systems  one  can  decrease  the  total 
system  cost  substantially  by  using  the  base-stock  policy  given  by 
(23)  rather  than  (24).  Finally,  we  note  that  the  critical  region 
is  larger  when  a  =  .8  than  when  a  =  .9. 

5.  CONCLUSIONS  AND  EXTENSIONS 

We  have  formulated  a  constrainted  NLP  to  solve  a  base-stock 
model  with  a  service-level  constraint  rather  than  shortage  costs. 
This  model  has  the  advantage  that  a  manager  can  more  easily  speci¬ 
fy  service  objectives  than  shortage  costs.  We  derived  first-order 
necessary  and  second-order  sufficient  conditions  for  an  optimal 
base-stock  policy.  For  the  special  case  of  zero  lead-time  and 
exponentially  distributed  demands,  we  derived  a  computationally 
efficient  algorithm  to  find  an  optimal  base-stock  policy.  We 
showed  that  for  some  inventory  systems  there  is  a  substantial  cost 
savings  when  using  this  algorithm  instead  of  the  independent  and 
identical  service  approach  that  is  often  used  in  practice. 

5 . 1  Unifo rmly  Distributed Demands 

One  can  do  analysis  similar  to  that  in  Section  4  for  the 
special  case  when  L^  *  0  and  demands  are  distributed  uniformly. 
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Given  that  the  demand  distribution  for  item  i  is  uniform  on 


a.  +  b. 

(a^.b^),  so  y.  *  j -  and  the 


2  ^b1  "  *1^ 
variance  is  oi  =  —  -yg-' — 


we  have  the  following  algorithm,  analogous  to  (23),  to  determine 
our  optimal  base-stock  policy  x. 


(1)  Solve 


0  for  X*  >  0. 


X*k. 


(2)  Set  v^  =  for  i  =  1,...,N. 


I  (28) 


i  i 


(3)  Set  xi  =  a..  +  /T2  o^a,  for  i  =  1,...,N.  J 


Analysis  similar  to  that  done  for  exponential  demands  could 
be  done  here.  We  note  that  although  the  uniform  density  is  a  two- 
parameter  density,  only  the  variance  parameter  is  used  to  deter¬ 
mine  the  v..  in  algorithm  (28).  The  reason  for  this  is  that  the 
mean  parameter  is  a  "shift  parameter,"  that  is,  its  only  effect  on 
the  density  or  cdf  is  to  move  it  right  or  left.  It  is  easily  seen 
from  (13)  that  a  shift  parameter  never  affects  the  value  of  v^ 
for  any  distribution.  Of  course,  the  value  of  is  shifted  as 
necessary. 


5.2  Possible  Directions  of  Research 

There  are  obvious  difficulties  involved  in  solving  the  equations 

★ 

for  the  first-order  necessary  condition  (13)  if  the  cdf  cannot 
be  written  in  closed  form.  A  direction  of  research  is. to  examine 
numerical  methods  for  solving  (13)  for,  say,  the  normal  or  gamna 
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distributions.  A  basic  difficulty  is  involved  in  approximating  the 
cdf  so  that  the  equations  (13)  can  be  solved,  and  so  that  the  solu¬ 
tion  obtained  is  a  good  approximation  of  the  true  policy  x.  ; 

According  to  Veinott  [2,  pp.  754-757],  the  generalized  base- 
stock  model  with  shortage  costs  still  has,  for  some  nonstationary 
parameters,  a  myopic  base-stock  policy  that  is  optimal.  A  simple 
model  of  this  sort  is  one  in  which  all  costs  are  stationary,  but 
the  demand  distribution  means  vary  in  time.  In  order  to  guarantee 
that  a  myopic  base-stock  policy  is  optimal,  it  is  necessary  that, 
for  each  item  in  the  system,  the  means  are  nondecreasing  in  time. 
This,  however,  is  not  strong  enough  to  guarantee  that  a  myopic  base- 
stock  policy  be  optimal  when  GS  service  constraints,  rather  than 
shortage  costs,  are  used.  Thus,  another  direction  of  research  is 
to  investigate  optimal  policies  for  this  model.  Since  they  will 
not  be  myopic,  this  seems  a  very  difficult  problem. 

A  more  promising  direction  is  to  introduce  a  set-up  cost  in 
the  model.  The  infinite  horizon  case  using  the  stationary  demand 
distributions  can  be  investigated.  This  model  seems  similar  to  the 
infinite  horizon  case  for  (10),  and  so  results  similar  to  those  in 
this  paper  may  be  forthcoming.  In  particular,  results  for  the 
exponential  demand  distributions  may  be  readily  derived. 


4 


J 
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APPENDIX 


Theorem  3:  For  fixed  a,  D  as  given  in  (27)  is  a  function 


y2  h? 

only  of  —  and 

u,  K7 


Proof:  We  note  that  for 


a,b,a,b  >0,  r  =  -  if  and  only  if 
D  b 


there  is  some  y  >  0  such  that  a  =  ya  and  b  =  yb.  Hence,  we 
must  show  that  if  =  yy^,  y^  =  yy2»  h.  =  6h^,  and  h^  =  6hp, 

then  D (y i  ,vi2 > ^ i  * ^2 ^  ”  D(y^  ,y2»h^  kt^).  Let  D  -  D(y-|  ,u2»h-j  ^2) 
and  6  =  D(y,  ^th.,^).  Also  let  and  be  given  by  (23) 
for  the  parameters  y.,h..  and  y^h..,  respectively.  Then  (27) 
implies 


0  =  1001 


2 

[a  +  log(l-a)]  ^  h.y. 

i  =  l  1  1 

2 

S  [v.-a  +  log(l-v.a)]h.y. 
i=l  1  ill 


and 


6  =  100  1 


2 

[a  +  log(l-a)]  Sh.-YV,- 
i=l  1  1 

"2 - 

£  [v.a  +  log(l-v.a)]5h.yy. 
i=l  1  ill 


since  y^  *  yy^  and  h^  =  6h..  Simplifying  yields 


28 


29 


and  so  since  X*  is  the  unique  solution  to  (29),  we  have 


it 

<5y 


=  x*  . 


Thus, 


r.  .  X*6yki 

Vi  a(x*6yki  +  "  Vi  * 


as  needed.  □ 
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